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The Elliott- Yafet (EY) mechanism is arguably the most promising candidate to explain the ul- 
trafast demagnetization dynamics in ferromagnetic transition metals on timescales on the order of 
100 fs. So far, only electron-phonon scattering has been analyzed as the scattering process needed 
to account for the demagnetization due to the EY mechanism. We show for the first time that 
the electron-electron scattering contribution to the EY mechanism has the potential to explain 
Ch , time-resolved magneto-optical Kerr efi'ect measurements on thin magnetic Co and Ni films, without 

^ ■ reference to a phononic "spin bath." 
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^ Current research in femtosecond magnetism is concerned with elucidating the fundamental mechanisms of light- 

O . induced spin dynamics as well as searching for potential applications in data processing [l, 2, 3]. Despite of important 
• experimental studies employing various time-dependent magneto-optical techniques, no consensus on a microscopic 
understanding of ultrafast magnetization dynamics in ferromagnets has emerged. Rather, demagnetization dynamics is 
typically described in the framework of the phenomenological three-temperature model. In this model, temperatures 
are assigned to the electron, lattice and spin "subsystems," and the exchange of energy (and spin) is driven by 
the temperature differences between the respective subsystems. Although the three-temperature model provides an 
intuitive picture of demagnetization, its relation to the microscopic dynamics behind the demagnetization is still an 
active field of research. 

The most popular candidate [3| for the microscopic process behind ultrafast demagnetization is a mechanism of 
^ ■ the EUiott- Yafet (EY) type [1]. In the EY mechanism, the demagnetization arises because, in the presence of the 
O spin-orbit interaction, spin is not a good quantum number, so that any momentum-dependent scattering mechanism 
changes the spin admixture when an electron is scattered from state to |fe -I- q). So far, the scattering process 
\ responsible for the EY mechanism has been assumed to be electron-phonon and electron-defect scattering in several 
—i • theoretical and experimental studies [1, H, 0, S, 0| ■ Unlike these papers, we apply, for the first time, an EY mechanism 
^ ' based exclusively on electron- electron Coulomb scattering to the ultrafast demagnetization in ferromagnetic metals. 
As a proof of principle, we demonstrate quantitative agreement for the demagnetization time and magnetization 
quenching between time-resolved magneto-optical Kerr effect measurements on Co and Ni, and the EY mechanism 
' based on electron-electron scattering. 

[ To resolve the electronic demagnetization dynamics on ultrafast timescales, we calculate the non-equilibrium, 
' momentum- resolved multi-band electron dynamics at the level of Boltzmann scattering integrals, and include the 
^\ ' carrier excitation process. We therefore do not include true electronic correlation effects beyond carrier scattering nor 
coherent effects due to the optical excitation process lol 11. 1^. 
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For the description of our general approach, let us assume that the electronic single-particle energies and wave 
functions fc), where the electronic band index fj, runs over majority- and minority-spin bands, and the vector 
rS momentum is labeled by k, are known from a band structure calculation. Then the Coulomb and dipole matrix 
1 elements can be calculated and used as an input for dynamical equations for the band and momentum resolved 
distribution functions, n^. From these, the total magnetization of the system is obtained by M = X]p=±fc ^M^-fc where 
= +1/2 for majority and — —1/2 for minority bands [2]| . The equation of motion determining the carrier 
distribution functions has the form [13| 
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For the electron-electron Coulomb scattering we use the Boltzmann equation in the form 
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where V is the dynamically screened Coulomb potential that depends on the initial and final states of the two scattering 
electrons fc) A; + q) and 1/^2, i + q) — > 1/^3, ^'^d hu = — e .!^ , . The optical excitation contribution in 

Eq. ID) is calculated by adiabatic elimination of the optical polarization [14| 
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where £J is the classical electromagnetic field, and the function g{hLu) (peaked around the central frequency of the 
excitation pulse) models the spread of photon energies which can induce electronic transitions via the dipole matrix 
element d^^(fc) between states A:) to Note that the electron-scattering based version of the EY mechanism 

is contained in Eqs. (Il|)~(l3]), if the Coulomb matrix elements include the spin-orbit interaction in the presence of the 
static lattice, so that scattering transitions change the average spin of the scattered electrons. The lattice effectively 
acts as a sink for the electronic angular momentum, which is "lost" from the electronic system by the spin non- 
conserving scattering processes described by Eq. The important Coulomb and dipole matrix elements can, in 
principle, be determined from ab initio treatments 0, and parameter-free results can be achieved by a dynamical 
solution of Eqs. (IT|)-(l3l). However, due to the numerical complexity of the fc-resolved Boltzmann scattering integral ([2]), 
we use a simplified model that contains parameters. We stress that the parameters and simplifications introduced 
below can be eliminated by using input from ab initio methods. We approximate the energy bands as spherically 
symmetric, = e|^|, and the screened Coulomb interaction as 

iX';,^(fc,£,9,c.)^/;v/^,^^ (4) 

where 
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The parameter a is roughly comparable to the a parameter introduced by Yafet fl^ and calculated recently for 
metallic fcrromagnets [7]. In Eq. ([¥]), v{q) is the bare Coulomb potential, and e~^(g,w) the dynamical inverse 
dielectric fmiction. For k, I, and q dependent factors /, this form can be shown to be valid if there are no short-range 
contributions to the Coulomb interaction [16], and we use this as an approximate explicit expression for the Coulomb 
matrix element of metals. Important dynamical screening effects are included via the Lindhard dielectric function 
e{q, uj). In semiconductors, it has recently been demonstrated that an approach closely related [l^ to the one presented 
he^e leads to a parameter free agreemeni for the spin dynamics in theory and experfment because quite accurate 
wavefunctions can be obtained using k ■ p theory. Finally, we assume that optical excitation connects only majority 
bands and minority bands with each other, respectively, and we approximate the strength of the optical dipole matrix 
elements by a momentum and band-independent constant d. Although this is a drastic oversimplification, especially 
in view of the hybridization between s and p bands, it is in the same spirit as the approximations introduced for the 
band structure and the Coulomb interaction: The dependence on the electron vector momentum should either be 
included in all these quantities, or modeled in a way that introduces the least amount of parameters in the model. In 
this paper, the electronic excitation after the optical pulse is therefore determined by the band structure, the central 
photon energy and width of pump pulse, as well as the flucnce. These quantities are used as input for the numerical 
calculations. 

Finally, we include the equilibration of the electronic system with the lattice via a relaxation time approximation 
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in Eq. ([T]). Here, denotes the Fermi-Dirac distribution of electrons in band /i at lattice temperature. Note that 
in our model the demagnetization solely occurs due to Coulomb scattering, and that the electron-phonon interac- 
tion only leads to thermal equilibration, which actually restores the ground-state magnetization. The equilibration 
times Tphon.Ni = 25 ps and Tphon.Co = 5 ps are extracted from experiment by fitting the remagnetization dynam- 
ics. This assumption is based on the experimental observation that the timescale for energy equilibration due to 
electron-phonon interaction, which includes the effect of heat diffusion in apragmatic way, is typically longer than the 
demagnetization time. 





FIG. 1: Band structure e|f of Ni (a) and Co (b) for majority (solid lines) and minority (dashed lines) electrons in F-X 
direction p^ . The Fermi energy Ef is set to zero. The arrows indicate a typical ultrafast demagnetization scenario: Electrons 
are excited by an ultrashort laser pulse (vertical light grey arrows), which does not change the total magnetization. They relax 
via intra- and inter-band scattering (grey arrows). The latter scattering process leads to depolarization of the electrons. 



On the experimental side, a variety of techniques are available to excite and detect electron-spin dynamics [l8l |. 
Here, we apply an all-optical strategy to trace the spin dynamics on femtosecond timescales. By means of the time- 
resolved magneto-optical Kerr effect (TR-MOKE) in the longitudinal configuration we excite the ferromagnet by an 
ultrafast optical pump pulse, and monitor the material response by a delayed and modified optical replica (probe 
pulse). The femtosecond pulses are generated by a Ti:Sapphire multipass amplifier with IkHz repetition rate. We 
use s-polarized 50 fs, 800 nm pump pulses at normal incidence, and s-polarized 50 fs, 400 nm probe pulses under 45°. 
The samples are thin polycrystalline ferromagnetic layers: a 15 nm cobalt film deposited on MgO by dc-sputtering, 
and a 15 nm Ni film deposited on Si by electron-beam evaporation. The Ni film is capped by a 3 nm Ti layer; another 
3 nm Ti layer acts as an adhesion promoter between the Ni film and the substrate. 

Qualitatively, the ultrafast demagnetization occurs in our model in the following way. The electronic distributions in 
the unexcited ferromagnet are assumed to be Fermi-Dirac distributions determined by the lattice temperature and the 
band structure. The majority and minority energy dispersions are spin split so that a non-zero magnetization exists 
in equilibrium. The ultrafast optical excitation process creates nonequilibrium electrons in bands accessible by the 
pump photon energy, and the electrons undergo intraband and interband Coulomb scattering processes. The driving 
force for the demagnetization are interband scattering processes between the optically excited electrons, which lead 
to the redistribution of electrons from majority to minority bands as long as the optically excited electrons are spin 
polarized. Remagnetization occurs due to equilibration at lattice temperature, because the groundstate magnetization 
is restored when the electrons settle down in the band minima. 

For the numerical calculations we use as input for a KKR-DFT-result [3| for the T-X direction, which is then 
used for the whole Brillouin zone as if the band structure were spherically symmetric. These dispersions are plotted in 
Fig[TJ For an experimental and theoretical study of Co and Ni band structures that also discusses the nomenclature 
of the bands, see Ref. The exciting laser pulse has a typical full- width at half-maximum of 50 fs, photon energy 
of 1.55 eV and the fluence is numerically adjusted to be in qualitative agreement with an estimate of the absorption 
and the observed magnetization quenching. The pump pulse excites electrons into initially empty states above the 
Fermi energy E-p, as modeled by Eq. ([3]). Some numerical results obtained from Eq. ([1]) for the time- and momentum- 
resolved electron occupation for Ni are shown in Fig. [2l Since the distribution functions contain all the information 
on the dynamics on the single-particle level, we use them, together with the band structure shown in Fig.[Tl to discuss 
the demagnetization scenario for Ni. Optical excitation by the ultrashort 1.55 eV, excitation pulse is only possible for 
transitions from the A J and bands to the A J (light grey arrows in Fig. [T]). Fig. ^a,) shows the nonequilibrium 
distributions created by the pump pulse. During and after the optical excitation of electrons in A\, electron-electron 
scattering processes redistribute the carriers in and between the bands. In Figs. [2] (c) and (d), the increasing number 
of electrons in Aj and Ag for positive time delays above k = 7nm^^ illustrates the dominant scattering pathways. 
Note that the scattering of electrons from majority to minority bands reduces the overall magnetization, as the 
electronic contribution to the expectation value of the spin is altered. The processes responsible for the ultrafast loss 
of magnetic order in Ni start more than 1 eV above the Fermi energy Ep) and mainly take place near the X-point. 
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FIG. 2: Dynamical distribution functions n'^ for four different bands. Electrons are excited from bands and Aj (a) into 
band AJ (b). Ultrafast demagnetization occurs by scattering into bands Aj (c) and Ag (d). (fe = 8.9nm~^x the index in the 
figures) 

Remarkably, the demagnetization is almost completely dominated by the two transitions mentioned above, with the 
nonequihbrium scattering dynamics taking place over more than 100 fs, cf. Figs. [2](c) and (d). Electrons that are not 
scattered out of band A| at high energies above the Fermi level accumulate in states close to the Fermi-level because 
energy- and momentum conservation requirements make out-scattering processes inefficient. For completeness, we 
mention that band Aj does not play an important role in the demagnetization dynamics of Ni. An analysis of the 
electronic occupation in the different bands for Co along the same lines leads to the scenario as depicted in Fig. [TJb). 

In Figure m we plot both the signal obtained from TR-MOKE measurements and the calculated signal for Ni and 
Co. The magnetization quenching for Ni is significantly stronger than for Co. In our model this is explained by the 
band structure in combination with the optical excitation process, which yields a more efficient carrier excitation with 
800 nm photons: in Ni mainly majority electrons are optically excited (see Fig. [ij so that all interband scattering 
processes above Ey lead to demagnetization. To obtain quantitative agreement between theory and experiment, we 
assume the same laser fluence for both materials, and use the Elliott- Yafet factor a, introduced in Eq. ([4]), as a single 
fit parameter. We obtain aco = 0.15 and am — 0.3 together with the demagnetization times of tqo = 215 fs and 
TNi = 200 fs. The general trend aco < c«Ni and the order of magnitude compare well with recent ab initio results for 
the a parameter [3]. The results in Ref. provide only a qualitative check for our fit parameters, because the ah initio 
results depend on the band-structure region, over which the wave-function coefficients are averaged. Last, but not 
least, we stress that the band structure properties influence the microscopic dynamics sufficiently strongly to make it 
impossible to fit the Co measurements in Fig. [3] using the Ni band-structure and, vice versa. 

In conclusion, we applied, an Elliott- Yafet mechanism based on electron-electron Coulomb scattering for the ex- 
planation of the optically induced ultrafast demagnetization in the ferromagnetic transition metals Ni and Co. The 
electronic demagnetization in our model occurs through the Coulomb interaction, which is spin-diagonal for free 
electrons, in the presence of the spin-orbit interaction. Modeling the optical excitation process and the scattering 
dynamics by Boltzmann scattering integrals for the momentum-dependent dynamical distributions functions in the 
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FIG. 3: Normalized TR-MOKE rotation (grey line) and calculated magnetization dynamics for Ni (blue line) and Co (red line). 
Assuming the same laser fluence, the choice of parameters qco ~ 0.15 and oni = 0.3 yields good agreement between theory 
and experiment. The demagnetization times are rco ~ 215 fs and tni ~ 200 fs. In the calculation, the stronger quenching of 
the magnetization for Ni is due to band structure effects and not very sensitive to the numerical value of a. 

various bands, we described a model that can make parameter-free predictions from input based ab initio band struc- 
ture results. In this paper, we made simplifying assumptions for the Coulomb and dipole matrix elements, which 
led to the introduction of two parameters, the EY-parameter a and the flucncc, which are fixed by comparison with 
experiment. Even though the origin of ferromagnetism is an interaction effect of many electrons, the good agreement 
between theory and experiment presents evidence that the ultrafast demagnetization dynamics in ferromagnets can 
be understood in a single-particle picture including electron-electron scattering and the spin-orbit interaction. Unlike 
other EY based models, we do not take into account the electron-phonon interaction for the demagnetization dynam- 
ics. Our goal was to present a proof-of-principle that ultrafast demagnetization can occur due to an EY mechanism 
based on electron-electron scattering. The question whether the electron-electron or electron-phonon interaction pre- 
dominates should be investigated further. Theoretically this can be achieved using the tools developed in this paper 
with the help of ab initio band structure calculations. Experimental strategies to separate the two contributions 
should be developed concurrently. 

We acknowledge support from the Graduiertenkolleg 792 "Nonlinear Optics and Ultrafast Processes" and SPP 1133, 
and from the NIC Jiilich through a CPU time grant. We thank Christoph Doring for sample preparation, and 
M. Fahnle (MPI Stuttgart) for helpful discussions. 
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